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A significant accumulation of matter in solid 4 He observed during the superflow events, dubbed as 
the giant isochoric compressibility (or the syringe effect), is discussed within the model of dislocations 
with superfluid core. It is shown that solid 4 He in a contact with superfluid reservoir can develop a 
bistability with respect to the syringe fraction, with the threshold for the bias by chemical potential 
determined by a typical free length of dislocations with superfluid core. The main implications of 
this effect are: hysteresis and strongly non-linear dynamical behavior leading to growth, proliferation 
and possibly exiting from a crystal of superclimbing dislocations. Three major channels for such 
dynamics are identified: i) injection and inflation of the prismatic loops from the boundary; ii) 
Bardeed-Herring generation of the loops in the bulk; iii) helical instability of the screw dislocations. 

It is argued that the current experiments are likely to be well in this regime. Several testable 
predictions for the time and the bias dependencies of the dynamics are suggested. 

PACS numbers: 67.80.bd, 67.80.bf 


I. INTRODUCTION 

The superflow through solid 4 He observed first in the 
UMASS groupP and then confirmed by other groupiPEI 
is now at the focus of the experimental and theoretical 
efforts in the field of superfluidity and quantum crystals. 
One of the striking features is the syringe effect (or the 
giant isochoric compressibility^. In its essence, a solid 
exhibits the response on external chemical potential ap¬ 
plied at a point, practically, the same way as liquid does 
- absorbs or expels a macroscopic fraction of atoms. 

As it has been suggested in RefP, this effect can be 
associated with the so called superclimb of edge disloca¬ 
tions - the climb supported by the superfluid transport 
along dislocation core. The unusual feature of this sce¬ 
nario is that the linearized isochoric compressibility of a 
solid permeated by a network of dislocations with super¬ 
fluid core is independent of density of the superclimbing 
dislocations and is, instead, determined by the dimen¬ 
sionless parameter - the asymmetry between lengths of 
superclimbing and non-superclimbing parts. This implies 
that the effect is strongly non-perturbative, that is, it 
cannot be treated as a small correction with respect to 
dislocation density. In particular, as shown in RefP, the 
linear isochoric compressibility of a symmetric network 
is essentially the same as that of a liquid. Here we show 
that, even if the initial network is strongly asymmetric in 
favor of the non-superclimbing (superfluid) dislocations, 
there are scenarios still leading to the giant isochoric com¬ 
pressibility. 

According to the suggestion^ the vycor ’’electrodes” 
are creating a contact between superfluid reservoirs and 
a preexisting static network of dislocations with super- 
fluid cores. Such dislocations are characterized by Burg¬ 
ers vector along the main symmetry axis and can be of 
two basic types - screv^ and edge^ (as well as of the 
mixed type). The edge part of the network can execute 
superclimb responsible for the syringe effect. An exam- 



FIG. 1. (Color online) A forest of screw dislocations con¬ 
taining edge superclimbing segments. Dashed and solid lines 
indicate pure screw and pure edge segments, respectively, all 
characterized by the Burgers vector along the hep axis. 


pie of such a network with combined segments is shown 
in Fig. |T| 

There is an alternative to the ’’preexisting static net¬ 
work” scenario - a dynamical network which is created 
and disrupted by the external bias p. Here the analysis 
of the superclimb is extended beyond the linear response 
considered in RefP, and it is shown that a segment of 
a rough superclimbing dislocation is unstable with re¬ 
spect to its unlimited growth, if the bias by chemical 
potential p exceeds a threshold p c which is inversely pro¬ 
portional to a length L of the segment. In high quality 
crystals typical value of L can be as large as several pin 
or a fraction of mm or even reach a sample size. Thus, 
the threshold can be macroscopically small, so that it 
may well be exceeded by several orders of magnitude in 
the experiments 1 3 . The reason for such a generic sit¬ 
uation is that a number of the conducting pathways is 
oc 1/A 2 , so that, in order to detect the superflow, the 
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FIG. 2. (Color online) A prismatic loop (solid line) gener¬ 
ated by an edge segment AB, from Fig. [I] according to the 
Bardeen-Herring mechanism: An originally straight edge seg¬ 
ment AB (solid horizontal line) bows under the bias (dashed- 
double-dotted line). Then, further bowing results in the over¬ 
hangs (dashed-dotted line). The points C,C’ in the overhangs 
approach each other and eventually the whole prismatic loop 
detaches from the points A,B. 


bias /i needs to be increased at least as /r oc L 2 which 
leads to fj, » /i c ex L~ l . 

It will also be shown that a straight screw dislocation 
(which cannot perform superclimb) with superfluid core 5 
can develop a helical instability under the bias so that 
the edge-like rim is formed and, accordingly, the syringe 
effect will also be induced. The threshold for this insta¬ 
bility has the same dependence oc 1 /L on length of the 
screw dislocation. [Helical screw dislocations have been 
first observed in silicon at high temperature^. 

The instability has two important consequences: First, 
rough superclimbing segments of dislocations pinned in¬ 
side solid 4 He bulk can generate prismatic loops upon 
bias /r in a manner very similar to the Frank-Reed source 
of gliding dislocation loops under shear stresW The dif¬ 
ference is that the generated loops during the superclimb¬ 
ing instability carry extra matter or vacancies. [In this 
case the instability should rather be called as Bardeen- 
Herring. A typical diameter of the loop is determined 
by the original length of the edge segment L = L 0 - This 
process is schematically shown in Fig. [2j Second, super¬ 
climbing dislocations existing at a solid-vycor boundary 
can proliferate into the bulk upon applying the bias, so 
that a percolating network of superfluid pathways is cre¬ 
ated even if it didn’t exist originally. This process is 
shown in Figs. |3|4| Both effects are symmetric with re¬ 
spect to the sign of the bias /x. In one case an additional 
matter is injected into the solid in the form of parts of 
extra basal planes, and in the second existing basal 
planes are being dissolved, that is, vacancies are being 
injected instead. 

The presented analysis is conducted at the level of a 
single dislocation. It ignores how pinning by 3 He impu¬ 
rities or crosslinks with other dislocations may affect the 
dynamics of the instability. It is clear that the insta¬ 
bility may also result in the dislocations exiting the solid 



FIG. 3. (Color online) Sketch of the growing superclimbing 
dislocation (dashed line). The area to its left indicates either 
building of extra plane of atoms or dissolving of the existing 
plane (not shown) between the upper and lower ones. 

from its edges. Several growing segments may also merge 
or recombine. These processes as well as the interaction 
between superclimbing and basal plane gliding disloca¬ 
tions are also not considered. In some sense the analysis 
presented here is limited by low density of superclimb¬ 
ing dislocations so that there is some reasonably long 
time during which the dynamics can be treated within 
an approximation of a single dislocation segment. This 
situation is different from the linearized approactf® where 
the main assumption was that a typical distance between 
superclimbing segments is of the same order as a typical 
length of the each segment. The single loop strongly 
non-linear dynamics considered here relies on a different 
limit - that is, a typical distance between superclimbing 
segments is much larger than Lq. 

The injection of the dislocations from the vycor-solid 
boundary and the Bardeen-Herring type loop generation 
as well as the helical instability of screw dislocations re¬ 
sult in the syringe effect. The dynamics of the instabili¬ 
ties, however, turn out to be different: while in the case 
of the boundary instability the injected dislocation can 
grow to sizes far exceeding its original length L = Lq 1 
in the case of the Bardeen-Herring type instability the 
generated loop radius R is of the order of L 0 ~ it is the 
number of loops that is changing. The helical screw dis¬ 
location can also generate loops in a manner similar to 
that proposed in RefP. [The detailed study of this effect 
in the context of the superfluid core will be conducted 
elsewhere]. 

II. GROWTH INSTABILITY OF ROUGH 
SUPERCLIMBING DISLOCATION 

Let’s consider one segment of a superclimbing disloca¬ 
tion of some initial length L = L 0 . Such a segment can 
be at a crystal boundary or be a part of the superfluid 
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FIG. 4. (Color online) Geometry of the growing superclimb¬ 
ing dislocation (dashed line) sketched in Fig. [ 3 ] The end¬ 
points of the dislocation are in a contact with the superfluid 
reservoir.The (solid) area under the curve indicates either the 
injected extra matter leading to the formation of new basal 
layer (actually two of them in hep ) or a removed part of the 
existing one. 


network. If the bias /x is applied, the segment will bow 
due to an extra matter delivered through its ends. Such 
bowing occurs in the basal plane while the Burgers vector 
b is perpendicular to the plane, that is, along the high 
symmetry axis. This process is schematically shown in 
Fig- HI as an arc protruding between two basal layers. 

At this point, let’s specify what the bias pi is. An in¬ 
crease of chemical potential of the superfluid reservoir pii 
either by applying pressure 1 3 or through the Fountain 
effeclP creates a difference pi = pi s —pii < 0 between chem¬ 
ical potentials of the solid pi s and the liquid. As a result, 
an additional matter can be injected into the solid in a 
form of growing pairs of basal plane layers. The bound¬ 
ary of one pair of layers is the superclimbing dislocation. 
If pi > 0, an existing pair of layers is being dissolved. Its 
boundary is also a superclimbing dislocation. The added 
or removed part of the planes is shown by a colored solid 
area under the arc in Fig|3] 

It is important to realize that imposing any finite value 
of pi (that is deviation from the equilibrium value) results, 
strictly speaking, in the instability. This can be under¬ 
stood from simple energy balance: the energy gain due 
to the bowing SEt = \pi5N\, with SN being a number of 
atoms delivered through the core to support the bowing, 
always exceeds the energy SE CI oc (T — T 0 ) due to the 
core length increase from T 0 to L for large enough L be¬ 
cause SN is given by the area swept during the bowing. 
Thus, for large enough L the energy gain due to the bow¬ 
ing always dominates the energy loss due to the length 
increase. However, in the limit pi 0 the (meta)stability 
is protected by a macroscopic energy barrier. This bar¬ 
rier vanishes if \pi\ exceeds a threshold ~ 1/To so that the 
absolute instability develops. As the estimates provided 
later show, the actual experiment^®] appear to be well 
in this regime. 


Let’s estimate the threshold value pi c for the bias. The 
energy of the dislocation per its unit length is given by 
shear modulus G and Burgers vector b as e c « Gb 2 / 47r, 
so that 5E CI « e c (T — L 0 ). The energy gain SEt scales 
by the area ~ \pi\(L/b) 2 swept by the bowing dislocation. 
Thus, equating one to the other gives pi c as 


Me 


Gb 4 

To 


(1) 


The value of b in Eq. 0 is the Burgers vector along the 
hep axis, that is, b = \J 8/3a where a is the interatomic 
distance a ~ 3.5A. As was found in the simulation^, this 
dislocation splits into two partials with 6 —>■ 6/2 and the 
fault in between. Thus, effectively, b is reduced by a fac¬ 
tor of 2 so that the actual threshold (|T]) becomes lower 
by a factor of about 2 4 . In the following discussions I will 
ignore this peculiarity of the structure, which can only 
modify the numerical coefficient without affecting the de¬ 
pendence pi c oc 1/T 0 . Accordingly, in all the following es¬ 
timates the value of b will be taken as b ss a ~ 3.5A and 
the core splitting will be ignored. 

The above relation can be supported by a more quan¬ 
titative analysis. In the quasi-static situation, so that pi 
is the same over the whole dislocation length, the bowing 
dislocation takes the shape of circular arc characterized 
by the base T 0 and the angle a as indicated in Fig. [4] 
At small \pi\ the center of the circle is outside the crystal 
so that the circle radius R >> Lq and ao —> 0. As the 
arc grows, ao eventually reaches 7 r and R decreases to 
R = To/2 and then ao 27 t so that the center of the 
circle enters the crystal and R starts growing to become 
R » To . The total energy of such a configuration can 
be represented as 


E = ^-R ■ a 0 - -yyl? 2 • (ao - sina 0 ) (2) 
47T b z 

where the radius R of the arc is determined by T 0 and 
the angle ao as To = 2sin(ao/2)I?. This expression indi¬ 
cates that the dislocation is absolutely unstable toward 
inflation R — > 00, ao —> 27r for arbitrary small \pi\ simply 
because the second term is negative and can dominate 
at large enough R. There is, however, an energy bar¬ 
rier to overcome before the instability develops unless \pi\ 
exceeds some critical value. Let’s consider a specific sit¬ 
uation when the end points of the dislocation are pinned 
by a contact with a superfluid reservoir so that To is 
fixed. Then, the energy ([ 2 ]) becomes a variable of the 
angle ao only: 


E(Lq, ao) — Eq 


a 0 sin(a 0 /2) - g ■ (a 0 - sina 0 ) 
sin 2 (a 0 /2) 


( 3 ) 


where E 0 « Gb 2 L 0 /8n, g = n\pi\L 0 /(b 4 G). As the analy¬ 
sis of the function ^ indicates, for g « g c = 0.5 there 
is a metastable minimum at a « Ag « 1 followed by 
a maximum at larger ao- At g = g c both,the minimum 
and the maximum, coincide at ao = 7r which is the inflec¬ 
tion point indicating ending of the metastability domain, 
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so that at g > g c the dislocation becomes unstable to¬ 
ward unlimited inflation ao —► 2n. As can be seen, the 
threshold g = 0.5, that is, 


Pc = 


Gb A 

2nL 0 


( 4 ) 


is consistent with the estimate 0 - 

In a general situation one should expect a distribu¬ 
tion of the lengths L$ so that some segments remain in 
a metastable equilibrium and some are overcritical. In 
what follows such a distribution will be ignored and it 
will be considered that there are M segments of some 
length L 0 in a solid affected by the bias /i. Typical values 
utilized in the flow experiment^ when the syringe effect 
was observed, are in the range p = 5 • 10 -4 — 5 • 10 _3 K 
(which corresponds to 0.001 — O.OlJ/g in units of Ref. 1 ). 
This implies that the lengths of the critical segment 
Lq > 1 — 5 pm (for a typical G ~ lOObar and b ~ 3.5A). 
The expected density of dislocations in high quality 4 He 
crystals is at the level of ~ 10 4 — 10 6 cm -2 as found in 
RefP, which implies that the actual lengths L 0 of free 
segments are about a factor of 10-100 longer than the 
above estimate. In other words, the experimental range 
of p used in Ref.® appears to be well above the threshold 
0 (or 0). At this point it should be mentioned that 
the dislocation density valued are more relevant to the 
glide effect than to the superclimb. Nevertheless, this 
order of magnitude estimate can be used as a figure of 
merit. 


A. Helical instability of the screw dislocation with 
superfluid core 


As found in ab initio simulation^, screw dislocation 
in solid 4 He with Burgers along the hep axis has a su¬ 
perfluid core. If this dislocation is straight, there are no 
edge-type segments on it, and, therefore, it cannot per¬ 
form the superclimb. Here it will be shown that, if such 
a dislocation is biased by chemical potential similarly to 
the edge segment discussed above, it will become unsta¬ 
ble toward forming a helix with its axis parallel to the 
original orientation of the dislocation. Such a helix has 
the edge-type rim and, thus, it can be a cause of the 
syringe effect. 

Let’s consider a screw dislocation of length L oriented 
along the z-axis (that is, the hep axis). Then, a position 
of the core can be described in the cylindrical coordinates 
by the radial distance r(z) from its original position r = 0 
(in units of b) as well as by the azimuthal angle 9{z). The 
energy consists of two terms: the work done pAN by the 
bias p to accumulate some amount of matter AN due to 
creating the edge-type rim and the energy ~ e c due to 
the core length increase : 

E s = £ dz \^Y~d z e + |p z rf + r 2 (<9 z 0) 2 ]| ,(5) 


where it was taken into account that the additional mat¬ 
ter per unit length of the core is dAN/dz = j s r 2 d z 9/2 , 
with 7 S = ±1 being the chirality (handedness) of the dis¬ 
location. In other words, the total amount AN is given 
by the projection of the helix on the basal plane times 
the number of the complete turns. The ’’sign” of the 
matter accumulation depends on : if the screw is right 
handed, y s = 1, and the helix is right handed, d z 8 > 0, 
the solid mass increases, that is, AN > 0. Similarly, 
AN > 0 for both the screw and the helix being left 
handed (y s = —1, d z 0 < 0). Conversely, the amount 
of the syringe matter becomes negative if the chiralities 
are opposite to each other. Eventually, it will be seen 
that the sign of the syringe fraction does not depend on 
the screw chirality and is solely determined by the sign 
of the bias p as AN ~ —p. 

As a specific choice of the boundary condition, let’s 
presume that this dislocation is pinned at its both ends, 
that is, r(0) = r{L) = 0. Then, the variation with respect 
to 9 gives the equation 

( 6 ) 


where the boundary condition is taken into account. Its 
substitution back to Eq. 0 results in the effective energy 
of the dislocation as 


E = 



P ) 2 

8e c 


d z r ) 



( 7 ) 


This expression features an instability toward unlimited 
growth of r. At small p, similarly to the case of the 
edge dislocation, the solution r = 0 is a metastable one. 
As Eq.Q indicates, there is a difference with the edge 
dislocation case - the screw does not show any linear re¬ 
sponse of bowing in the limit p —> 0. There is, however, 
a threshold p s such that at \p\ > p s the absolute insta¬ 
bility toward r —>• oo develops. In order to find how p s 
depends on L and e c it is enough to perform elementary 
estimates: r << L changes on the scale of L so that the 
total elastic energy is ~ e c r 2 /L. As long as the bias en¬ 
ergy ~ p 2 Lr 2 /e c becomes of the same order, the solution 
r = 0 becomes absolutely unstable. Thus, 


Ps ~ 



Gb 4 
~L' 


( 8 ) 


which is essentially the same condition as (jT|) . At 
the threshold the helix is described by the total angle 
|a| ~ p s L/e c ~ 1, and at \p\ » p s this angle becomes 
\p\/p s » 1. 


B. Collective elastic effect in a bulk network of 
superfluid dislocations 

The injection of extra matter (or vacancies) under the 
bias p is limited by the compression elastic modulus K of 
a sample, so that the system stabilizes at some finite den¬ 
sity of extra matter delivered through superclimb. The 
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above estimate |T]) (or ©) obtained for a single dislo¬ 
cation does not take into account this effect and implies 
that an inflating loop can reach a sample size. In reality, 
the generation must stop after the overall density change 
compensates for the bias y. 

Let’s presume that the extra matter resides in M 
dislocation segments of length Lq which bowed by an 
amount y each. Such bowing results in the extra mat¬ 
ter AN ~ LyM/b 2 added to (or subtracted from) the 
solid. There is the corresponding compression energy in¬ 
crease E e « K(AN/N) 2 fl/2, where N stands for the 
total number of particles in the bulk of the volume f 1 
affected by the injection. The chemical potential energy 
gain and the energy loss due to the core length increase 
are ~ yLyM/b 2 and ~ Gb 2 y 2 /L, respectively. Thus, the 
total energy change as a function of y, M, L becomes 


E 


LMy 

b 2 


y + 


'Gb 2 M 

L 


I\L 2 b 2 M 2 ' 

n 



(9) 


where the dimensionless numerical coefficients are omit¬ 
ted.If \y\ « L, the value of L can be set to L ~ Lq and 
the minimization in y gives 


y 


LU(Gb*) 

1 + KL%M/(QG) 


( 10 ) 


As mentioned above, this solution is actually a 
metastable one, which, however, is protected by expo¬ 
nentially long waiting time in the limit y —> 0. The 
bowing determined by (101 corresponds to the syringe 
fraction AN/N ~ MLybJXl which for the case depicted 
in Fig. [l] where M ~ Ll/(L 2 L Z ) can be written as 


AN ~ L 3 /(L 2 L Z ) _y_ 

N ~1 + KL 3 /(GL 2 L Z ) Gb 3 ' 


If L s ~ L z ~ L, that is, for a uniform network of the 
superclimbing dislocations, this fraction becomes 


AN ^ y 
~N~ ~ (K + G)b 3 ’ 


( 12 ) 


which constitutes the giant isochoric compressibility. In 
the limit Lq << L 2 L Z , the syringe fraction becomes re¬ 
duced within the linearized approach as 


AN ^ Lq y ^ ^ y 
N ~ L 2 s L z Gtf << Gb a’ 


(13) 


If the bias y exceeds the threshold 0, the bowing of 
the edge segments cannot be treated in the linear approx¬ 
imation anymore. In order to find the syringe fraction in a 
generic situation one can use Eq. 0 where the substitute 
y ~ L is made and, instead of L, the fraction Ni ~ L 2 /b 2 
generated by one segment is used as a variable. Then, 
Eq.@ becomes: 

E « — \y\MNi + G&My/Th + N 2 , (14) 


where N\ is taken as a positive value featuring either 
extra matter delivered to (y < 0) or taken out from (y > 
0) the solid. 


For small enough M this function of N\ features a 
maximum at Ni = N mx ~ (Gb 3 /\y\) 2 and then a stable 
minimum at = N eq where 


Neq 


MKb 6 ’ 


(15) 


This minimum corresponds to the syringe fraction 
AN/N = N\Mb 3 /fi, that is, 


AN _ y 
~N~ ~ Kb 3 ' 1 


(16) 


where the condition N mx « N eq , that is, 


\/A > ^ 



(17) 


must hold. In the case of M bulk segments distributed 
uniformly in. e.g., the situation depicted in Fig. [T] the 

condition |l7| becomes \y\ > (G 2 K/L 2 L z ) 1 ^ b 4 . Thus, 
if Lq is the smallest length scale and \y\ obeys (l]), the 
system is guaranteed to be unstable, with the equilibrium 
(15) to be determined by the bulk elastic energy E e . The 


fraction (15) corresponds to the limit |/x| >> yb- 

In the case of the bulk structure shown in Fig. [l] this 
fraction Ni is residing in several prismatic loops Ni p of 
a radius R ~ L 0 generated by each edge segment. This 
number can be estimated as N eq b 2 /Lq, that is, 


N, 


i P 


\iA L l L z 

Kb 3 L 2 


> 


GLjL z 

“ALf 


2/3 


» 1, 


(18) 


where the condition <0 as well as that L = Lq is the 
smallest distance among L s , L Z: L in Fig.[l]are taken into 
account. 

Thus, in the overcritical regime the equilibrium frac¬ 
tion (16) is always of the same order as in the liquid 


- even if the linearized response predicts much smaller 
values |l3|). 


C. Collective elastic effects due to the boundary 
instability 


In the case of the vycor-solid boundary, the analysis 
should be performed separately because seeds of the un¬ 
stable dislocations are residing at the boundary. Accord¬ 
ingly, M in this case is rather a surface than the bulk 
quantity. Then, in the estimate of the bulk deformation 
energy the affected volume becomes Ll ~ LS, where S 
stands for the area of the vycor-solid boundary. Accord¬ 
ingly, the extra fraction of the injected matter/vacancies 
is A A” ~ L 2 M/b 2 and N ~ f l/b 3 so that 


AN LMb 

1st 7T’ 

and the elastic energy takes the form 
„ KM 2 b 2 L 3 

Ep_ ~ 


2 S 


(19) 


( 20 ) 
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This dependence ex L 3 should be contrasted with the 
elastic term oc TV 2 ex L 4 in the case of the bulk instability 
as represented in Eq.(14). The total energy in terms 
of Ni (that is, the extra matter due to one segment) 
becomes 

E sa — \n\MNi + Gb 3 M^/Ni + ^ 

This form as a function of Ni (or the segment length L) 
can also have two extrema - a maximum followed by the 
minimum as Ni grows. It happens when 


> Ms 


/ GKM\ 1/2 4 


V s 


- 


b . 


( 22 ) 


This condition for the surface bistability should be com¬ 
pared with the bulk one ( fl7| ). The value of M/S is de¬ 
termined by typical distances between the boundary seg¬ 
ments along the basal plane and along the hep axis 
r z as M/S ~ l/r\,r z . Thus, if Lq < y/ri,r z , the con¬ 
dition (22) is guaranteed to be satisfied as long as the 
instability condition ([!]) holds. If |/i| >> y s , the equi¬ 
librium is determined by the first and the last terms in 
the energy ©. It corresponds to the typical equilibrium 
length (obtained from the minimization of E in (21) with 
respect to N\ ~ L 2 ) as 


L = L e 


ImI nr z 
Kb 4 ' 


(23) 


It is interesting to note that this length becomes of the 
order of a sample size (~ 1 cm) for the smallest values 
of the bias ~ 5 • 10~ 4 K used in Ref.® (or y ~ 10 _3 J/g 
in units of Ref. 3) , if Vb,r z are of the order of the vycor 
diameter ~ 1mm. This value, however, drops quite fast 
with the product ri,r z « 1mm 2 . In other words, if there 
are only few seeds of superclimbing dislocations at the 
solid-vycor boundary, the instability guarantees that the 
new pathways will reach the other electrode. Conversely, 
if there are many such seeds, the elastic energy increase 
due to the injection will stop the syringe effect close to 
the boundary. 

Generically, a system of superclimbing dislocations can 
feature two minima with respect to the syringe fraction. 
These minima are separated by a barrier, and, therefore, 
the hysteresis phenomenon should be anticipated with 
respect to the bias as long as the condition @ (or (HO 
is satisfied. 


D. Renormalization of the chemical potential 

The instability of a single superclimbing dislocation is 
eventually stabilized by the increase of the bulk elastic 
energy due to finite density of the injected matter (or 
vacancies). This corresponds to the renormalization of 
the difference of chemical potentials y between solid and 
liquid from its initial value to zero. This renormalized 
value p can be obtained from the expressions of the total 


energy, Eqs.( 14|2l |, as p = d{E/M)/dN\. Keeping in 
mind the symmetry y —» — y let’s consider y < 0, that is, 
that the potential of the liquid in vycor is higher than in 
the solid so that extra atoms enter the solid. In the case 
of the bulk segments (as in Fig. |T|) the differentiation of 
the energy (14) results in 


/« = ),(»,) = P i(24) 


and in the case of the vycor-solid boundary Eq.(21) gives 


~ _ /Ar x Gb 3 , r -i/2 3Kb 5 M i /2 

M = MW) = M + -j-N. 1/2 - ^ 1/2 


i ( 25 ) 


There are two roots of p = 0. At small enough M ( as 
presented in Eqs.( 17|22 ) and guaranteed by the generic 
condition (|T]) ) the first one N± « (Gb 3 /(2y)) 2 corre¬ 
sponds to unstable equilibrium, and the second one de¬ 
scribes stable one. It should be mentioned that the equi¬ 
librium characterized by small bowings (that is, Ni —>• 0) 
is not captured by Eqs.( 24|2~5 | written for already large 
bowings y ~ Lq. Thus, for all practical purposes this 
minimum can be viewed as corresponding to the energy 
E = 0 reached at Ni = 0. 

As discussed above, the practical values of y corre¬ 
sponds to the situation when N\ is to evolve from the 
unstable toward the stable equilibrium. This what is 
called above as the overcritical regime \y\ > y c , Eq. 0 
As will be seen, the overcritical dynamics exhibit strongly 
non-linear features before Ni approaches the vicinity of 
the stable equilibrium. 


E. Liquid-gas type transition and hysteresis 


It is useful to look on Eqs.( 14|2l 1 from a different per¬ 
spective. Small bowing of dislocations corresponds to 
Ni 0, that is to zero energy E = 0. There can exist 
another equilibrium solution dE/dNi = 0 characterized 
by finite N\. Thus, there is a value of chemical potential 
|/i | = yj at which two phases Ni « 0 and Ni finite have 
the same energies. This can be interpreted as a point 
of first order phase transition. For the case of the bulk 


system,Eq.(14) 


yi = 1.5 ( G 2 K^\ ' b 4 . 


(26) 


At smaller values of \y\ the second solution for be¬ 
comes metastable and at /i = y sp , where 

y sp = 2 _1 / 3 /i/ » 0.794/ij, (27) 

it vanishes. Thus, \y\ = y sp corresponds to the spinodal, 
that is, to the point where the bistability (and hysteresis) 
vanishes. 

Similar situation occurs for the case of the interface, 


Eq.(21). The transition occurs at \y\ = yj s , where 

,M\ 1/2 


yi: 


-*{ aK T) b ‘ 


(28) 
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and the hysteresis vanishes at \p\ = p sps , where 


73 

f-^sps - 2 


0.866/i/ s 


(29) 


The transition is not characterized by any underly¬ 
ing symmetry, and, to some extent, resembles 1st or¬ 
der liquid-gas transition, where density exhibits a jump. 
However, there is also a significant difference. The en¬ 
ergies contain essentially a non-analytical term 

~ 7iVi determined by the geometrical nature of disloca¬ 
tions. This term is always dominant at small Ni and is 
the reason for the energy barrier. Thus, in contrast to 
the standard liquid-gas transition, the syringe effect does 
not have a critical point where the first order transition 
ends. 


III. DYNAMICS OF THE INSTABILITY AT 
THE VYCOR-SOLID BOUNDARY 

Now let’s consider the dynamical aspects of the evolu¬ 
tion of the syringe fraction. In this work, the focus is on 
the dynamics of the edge segments. The dynamics of the 
helical instability will be analyzed elsewhere. 

A single loop dynamics during the instability stage is 
characterized by a short ballistic period during which 
phase slips and dissipation can be ignored. Then, as 
the superfluid velocity along cores reaches some terminal 
value, phase slips induce dissipation which is strongly 
non-linear in the velocity. Finally, once length of a grow¬ 
ing segment approaches equilibrium value (determined 
by the largest root of p = 0 in Eq.p5|), the linearized 
dynamics sets in. The analysis is conducted within the 
assumption that distance between dislocations is large 
enough so that there is enough time for a segment to 
grow to a length L which is much larger than the initial 
length Lq. Within this approach other growing loops are 
taken into account self-consistently through the renor¬ 
malization of the chemical potential given by the last 
term in Eqs.( 24|25 |. Practically, this means distances 
~ 10 -3 — 10 -2 cmfor dislocation densities 10 6 — 10 4 cm -2 . 


A. Ballistic growth from the boundary 

The analysis will be conducted for an almost circular 
loop of radius R (that is, ao ~ 27r in Fig. 4j. In general, 
the parametrization of the loop should include deviations 
of its shape from circle. This can be achieved by intro¬ 
ducing R as a function of polar angle a (as defined in 
Fig. |4| and time t. The Lagrangian can be written as 


C = 


da 


— <j)aR 2 


-£r 


(30) 


where the first term is the Berry contribution in units 
Tl = 1,6 = 1 , and a = ±1 is to be chosen depending 
on whether the matter is injected to [a = +1) or from 


(a = — 1 ) the solid; the term E\ is the energy of one loop 
which takes into account the elastic energy, that is, the 
last term in Eq. (21), attributed to one loop. In addition, 

the kinetic energy of the flow along the core P 2 dijdl > 
where 4> is the superfluid phase and p s stands for the 
superfluid stiffness and ^ R 2 + (d a R) 2 , should be 
taken into account. Thus, 


r 2ir 


E-i = 


da 


p s {d a 4>f 

2 dl/da 


i dl 

-MiVr + ec — 


kA^I) 


where k « Kb 5 M/S , e c = G6 3 /47 t and Ni as the 
amount of extra matter absorbed by one loop is N\ = 
/ Q 27r daR 2 / 2. As discussed above, the system is invariant 
with respect p —> —p. Thus, without loss of generality 
the value of p can be taken negative and a = +1 so that 
N\ describes extra matter added to the solid. There is 
also the boundary condition indicating that the disloca¬ 
tion is in a contact with the superfluid reservoir at its 
two endpoints, that is, <p(a = 0) = <f>(a — 2n) = <j)R, 
where 4>r is the phase of the reservoir, which can be set 
to zero. 

Within the simplified approach deviations from the cir¬ 
cular shape can be ignored, that is, d a R = 0. Then, a 
certain minimal assumption must be made about the spa¬ 
tial dependence of the phase <j> along the dislocation core. 
Given the symmetry of the problem, Fig. [4j and because 
the total current through the loop must be zero in the 
syringe regime, the current along the core ~ d(f>/da must 
be antisymmetric with respect to a = tt .Thus, the lowest 
non-trivial angular harmonic satisfying this requirement 
as well as the boundary condition is 


(f) = 4>o(t) sin . 


(32) 


Then, a substitution of this ansatz into Eqs.( 30|3l l, after 
performing explicit integration and variation in 0 q, gives 


(33) 


dR 2 _ irp s <t> 0 
dt 8 R 

This equation is essentially the statement of the continu¬ 
ity: the flux of matter through two ends of the growing 
loop controls the rate of the loop area change. The vari¬ 
ation of the action in R gives 


d<j> o 

dt 


= P, 


where 


p = 


7/i| 

2 


+ 


np s (/)\ 


ne c 


- 3t t 3/2 kR. 


(34) 


(35) 


32 R 3 2 R 

For \p | exceeding the threshold 0 and when R is yet far 
from the equilibrium, the dominant time dependence is 
determined by the first two terms in p. Then, the so¬ 
lution of Eqs.( 33|34|35 1 can be looked for in the form 
(f> o = At , R = Bt v with some unknown parameters 
A,B,v. This gives 


<h{t) = 


,; i ( ( ) = (7gM) 1 'V (36) 



















These dependencies describe the balistic stage of the sy¬ 
ringe effect far from the equilibrium. The accumulated 
fraction is given by Eq.|l9|) where the role of L is played 


by R from Eq. (36), that is, 


oc (p s |p|) 1 / 3 t 2/3 , 


(37) 


which corresponds to the local pressure variation ~ 
K\AN\/N. These fractional powers are specific to the 
dislocation superclimb and, thus, their experimental ob¬ 
servation would be a ’’smoking gun” for the effect. The 
question, though, is how long this ballistic stage can last. 

During the ballistic stage the flow velocity grows in 
time until some terminal velocity Vr is reached (at the 
dislocation end points). Then, frequent phase slips take 
place which convert kinetic energy of superflow into ex¬ 
citations. Thus, as an order of magnitude estimate, a 
typical time for the phase slip t ps can be taken from the 
ballistic stage — how long it takes to accelerate the flow 
from zero to, say, Vr ~ lOOm/s. The velocity profile 
along the growing loop is determined by the phase (32) 
as V(t,a) = d a (f>/(mR), that is, 


v,t ’ a> = 2^m cos,a/2> - 


(38) 


where m is 4 He atomic mass. Using the solution for <j >o 
and R(t) in this equation, one finds 

h / t \ 1//3 

V(t,a) = Vq cos(a/2), U 0 = — f — J . (39) 


in the standard units, where the time scale Tb is given by 


37 TfVb p 2 mb 

~^p7~ ~ ' KF’ 


(40) 


and n's d ' > ~ mp s is the superfluid linear density along 
the dislocation core. Using its value ni ld ^ ~ lA _1 found 
in simulations of screw dislocation 5 and b ~ 3.5A(so that 
%/mb ~ 50m/s), the terminal speed Vr (at a = 0, 27t) is 
reached at time tr ~ Tb, which for the lowest p value used 
in Ref . 1 gives ~ 1 ms. At this moment it looks unlikely 
that the available time resolution allows observing the 
time dependence during this stage. There is, however, 
an initial portion of the syringe effe ct w hich is on the 
experimental time scale (of minutesj s accumulated 
in a jump-like manner right after the bias p is imposed 
(or removed). This quantity will be discussed later. 

Close to the stable equilibrium point R —>• R eq (deter¬ 
mined by (j) = 0 ), where 


/? = H 

eq 6 ^/ttk ’ 

which is essentially Eq. ([23]) , Eqs.( 
earized R = R eq +£ with |£| << R eq . 
dynamics is oscillatory: 


(41) 


33|34| ) can be 


lin- 


The corresponding 


£ + — 0 ) — 


277 T 7 / 2 p s K 3 
4 p 2 


(42) 


Thus, scales as oc l/|/x|. These oscillations, how¬ 
ever, can take place only if their amplitude is small so 
that the velocity of the flow along the core remains much 
smaller than the terminal velocity. Thus, detecting them 
presents a significant challenge. It should also be men¬ 
tioned that the phase slips in the ohmic regime at fi¬ 
nite temperature T (see below) may make the oscillations 
overdamped. 


B. Dissiptaive stage 


At the experimental time scale of minutei li2 3 the dy¬ 
namics is dominated by diffusive (dissipative) processes. 
The nature of these processes is not exactly known. One 
possible scenario is that quantum phase slips assisted by 
thermal processes in the superflow along dislocation cores 
are responsible for the dissipation. The dependence of 
the flow rate vs bia^S is consistent with the picture of the 
phase slips in Luttinger liquid containing a weak link 10 11 . 
At the sam e ti me, the origin of a significant temperature 
dependence 1 2 - observed in the non-linear regime is not 
fully understood. In this situation, a phenomenological 
approach should be used in order to describe the loop in¬ 
flation. Specifically, in the ballistic dynamical equation 
(34) the l.h.s. can be rewritten in terms of the veloc¬ 
ity amplitude Vq in Eq.(38) as </>o —I mRV V Then, the 


effective friction rate j(T, IV)Do, with some friction coef¬ 
ficient 7 depending on Vq , T, should be added to the flow 
acceleration dV a /dt in Eq.(34). This transforms Eq.(34) 
into 


p ,; 


dV ° 1 \r ~ ^ 

hF +lV °~^R 

■n\p\ irm 2 p s Vo ne, 


- ^ - 3tt 3 / 2 kR, 
32 R 2 R 


(43) 

(44) 


where the Bernoulli pressure is now expressed in terms 
of the velocity rather than the phase. 

The growth rate of the loop area ~ R 2 is determined 
by Vo- The actual relation (stemming from the continuity 
equation) is exactly the same as in Eq. (33) where, how¬ 
ever, the phase is now expressed in terms of the velocity 

w 


ndR 2 

dt 


— 2p s IV, 


(45) 


where it is understood that V = V(a = 2ir) = — V(a = 
0 ) so that the matter is delivered symmetrically from the 
both ends of the growing loop of radius R» L 0 . 

Let’s now make a choi ce for 7 . According to the quan¬ 
tum phase slip scenariotP - 0 U energy of ID superflow is 
converted into excitations of the Luttinger liquid. Gener- 
ically, the quantum effects assisted by thermal excitations 
are characterized by power law dependencies of the phase 
slip rate t~^ oc with some C, > 0 determined by the 
Luttinger liquid parameter g. At zero T the flow velocity 
V controls the dissipation. Within the weak link situa¬ 
tion a phase jump by ~ 7 r occurs at microscopic distances 
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across the link which is reasonable to take as ~ b. This 
jump is accompanied by energy transfer ~ V between 
the link and the Luttinger liquid. Thus, Vo plays the role 
of temperature so that the dependence on Vo should be 
characterized by the same exponent t~Y oc |Vop, with 
the crossover taking place at some Vo = V T ~ Tb. Thus, 
as a single equation the friction rate in Eq. ([43]) can be 
represented as 


' ps 


'yVo = 7 0 Im[ 6 T + *V r 0 ] c , 


(46) 


with some coefficient 70. In the weak-link scenario 70 
is determined by frequent phase slips occurring at the 
location of the link. Thus, the associated time constant 
oc can be much shorter than a typical time-scale 
set by the period of Debye frequency in solid 4 He. This 
may essentially eliminate the ballistic stage for practical 
durations of the experiments. So, below the estimate for 
the loop inflation will be obtained under this assumption, 
that is, that the ballistic stage is too short to produce any 
significant syringe effect. 

The power ( can be empirically related to the power p 
observed in the flow rate vs bias dependence |Vj oc |p| p 
in Ref. 1 ! as 7V0 ~ |p| p so that 


C =P~ 


(47) 


According to Refs® ( = 2/g ~ 1 and the self-consistent 
result® gives ( = 2g — 1 . In what follows the power 
p will be used as a quantity measured directly in the 
experiments This power was found to vary in the range 
0.25 <p< 0.5. 


For \ V\ « Vr = Tb Eq.(46) implies ohmic regime 


7^0 = 7oIm[6T + iV 0 ] 1/p -> p-^oW^Vo. ( 48 ) 


Eqs.( 46|48] will be used below in the analysis of the 
loop dynamics in the long-time limit where the inertial 
part in Eq.(43) can be omitted. Then, far from the equi¬ 


librium at low T the dynamics is dominated by the first 
term in the brackets of Eq. (44). Thus, 


+n 


, -V 

\m'yoRJ 




in the non-linear regime (46), and 


„,p(bT)i-P~ 

F ° ~ m 10 R R 


(49) 


(50) 


in the ohmic regime (481. 


These expressions must be used in Eq.(45). Accord¬ 


ingly, in the non-linear regime the loop radius obeys 


dR 2 

dt 


2 p. 


\m'y 0 Rj 


(51) 


which implies R oc |p/7o| p A 2+p !(fp s ) 1 A 2+p ), or for the 
syringe fraction 

^ociiloc (My +P (/M)5 fe. ( 52 ) 


As Eq.(49) indicates, the flow speed actually drops with 
time as V oc t -1 /( 2 +p). Thus, eventually, the non-linear 
regime must change to the ohmic one characterized by 


dR 2 
dt 


2 P. 


pT 


i-p 


TO70 


R ’ 


which gives 


AN 

~N~ 


oc R oc T 3 


P s \p\t 

7o 


(53) 


(54) 


Eqs.( |52|54| are obtained under the assumption that 
the system is far from the equilibrium. If, however, it 
approaches the equilibrium, the last term in the brackets 


of Eq. (44) becomes important (with the second and the 


third ones still being irrelevant). This term stabilizes 
the system at the equilibrium radius R eqi Eq.(411. Close 


to the equilibrium the dynamics becomes linear in the 
deviation \R eq — i?| << R eq . As mentioned above, the 
time dependence would become either dissipative at high 
T or oscillatory as in Eq.(42). 


At this point it should be mentioned that the stabi¬ 
lization of the instability may also happen due to the 
dynamical rather than due to static equilibrium. Specif¬ 
ically, if R eq exceeds a system size, the stabilization is 
to be achieved by the balance of growing new loops and 
the loops exiting the sample. This picture essentially 
depends on sample geometry and size and will not be 
discussed here. 


C. The jump in the syringe fraction due to the 
ballistic inflation 


As discussed above, the ballistic stage may lead to a 
jump of the accumulated syringe fraction right after the 
bias is applied (or removed). Let’s estimate this fraction, 
first, for T = 0. In the dynamical equation ( [43] ) the dis¬ 
sipative part can be ignored as long as \V\ » 

Using the ballistic solution (39) in this estimate, the lim¬ 
iting time becomes 


1-P 
. 1 + 2 V 


Ps %Pp \d\ 


tbal 


7o 


3 p 

1+2 v 


OC 


2 ( 1 —p) 

l + 2p 


7o 


3p 

l + 2p 


(55) 


where the definition (40) of is used. A substitution of 


it into the ballistic fraction, Eq. (52), gives the jump as 


|AiV| 

N 


7o 


2p 

l + 2p 


1 — 2p 

\g\l+2p 


I Ml > Pc- 


(56) 


The value of p was found in Ref.® to be below 0.5. Thus, 
the ballistic jump is a decreasing function of the bias, 
provided it exceeds the threshold for the instability and 
the jump itself does not exceed the equilibrium syringe 
fraction oc |p|. [In this case, the last term in p, Eq.(44), 
should be taken into account which will change the bal¬ 
listic solution (39)]. However, as mentioned earlier, the 
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friction ” amplitude” 70 in the denominator may actually 
suppress the jump below the experimental resolution. 

Let’s now consider finite T. Comparing the accelera¬ 
tion rate with the thermal phase slips in Eq.(43) the bal¬ 


listic evolution takes place (before it is interrupted by the 
ohmic regime) as long as t is shorter than the smallest of 
either the ohmic dissipation time 7 c j“ 1 T ,1_1 / p or the time 
when the terminal velocity V = bT is reached. Clearly, 
at very s mal l T and p < 1 the latest dominates, which 
from Eq. (39) follows as t ~ t F oc 17 T 3 . Then, at longer 


times the evolution becomes essentially the same as that 
at T = 0 and leads to the estimates ( 55|56 ). However, at 
the experimental values of T and large , the estimate 
t < tr » 7%(7o 1 r 1 " 1 /*’) 3 is more appropriate for the time 
limiting the ballistic evolution. Then, a substitution of 
t F into Eq.(371 gives the jump as 


|AAT| psT^-P 

a 7o H ' 


(57) 


This dependence should be considered in the limit 70 —► 
00 , that is, that the maximum typical time for the phase 
slips ~ y ^ 1 is below (T/T 0 ) 2+ ^ p T b , where T 0 ~ IK is 
a typical temperature corresponding to the velocities ~ 
lOOrn/s. At this point, it is worth mentioning that the 
jumps of the syringe fraction have been observed in Ref. 2 . 
To what extent these can be interpreted in terms of the 
ballistic stage remains to be seen. 


the equilibrium p = 0 is reached. Thus, the accumu¬ 
lated fraction (far from the equilibrium) can be written 
as A N/N oc t/t FR » 1. Thus, the syringe rate dAN/dt 
becomes 


dAN 1 p s \p\ p 
dt “ t ffi K y 0 l 2 0 +p 

in the non-linear regime and 

dAN 1 p s T 1 - p_ Vl 

— , — oc - 00 - ^ - 

dt t FR IqLq 


(58) 


(59) 


in the ohmic one. After the bulk accumulated fraction 
approa ches the equilibrium value (16) the constant rates 
( |58|60 ) transform into the exponential diffusive slowing 
down. In contrast to the boundary instability in this case, 
the Bardeen-Herring mechanism is inherently dissipative 
and no oscillations are to be anticipated. 

For short loops the Bardeen-Herring cycle can occur 
in the ballistic regi me. In this case the time t FR can be 
estimated from Eq.(36) as t FR oc L^ 2 {p s \p\)~ 1 / 2 . Thus 
the rate becomes 


dAN 1 (p s \p\) 1/2 

dt t FR L y 2 


(60) 


V. DISCUSSION 


IV. THE BARDEEN-HERRING TYPE 
INSTABILITY 


While a dislocation injected from crystal edge to the 
bulk can grow up to R eq which is much larger than its 
initial size ( or even as large as sample size), a finite 
superclimbing segment inside a solid, e.g., in the case 
shown in Fig. [l] can generate loops only of a size of the 
order of its initial length. According to the Bardeen- 
Herring mechanisnP originally considered for gliding dis¬ 
locations and known as Frank-Reed instability®, an ini¬ 
tially straight segment bows under the bias, and eventu¬ 
ally the overhangs are created, Fig. [2] These overhangs 
merge together so that a circular (prismatic) loop of a ra¬ 
dius Rli which is of the order of initial length Lq of the 
straight segment, is created. This process is cyclic and 
is characterized by time t FR needed for the loop to grow 
until the overhangs (C,C’ in Fig. [ 2 ]) merge together so 
that the loop becomes separated from the main network. 
At this point what happens to this loop is not important 
it can,e.g., diffuse away or merge with newly created 
loops. 


An estimate for this time can be obtained from Eq.( 51) 


in the non-linear regime or from Eq. (53) in the ohmic 


regime, where the time t FR is found as a function of the 
loop radius R reaching the length of the order of the orig¬ 
inal segment length Lq. It is natural to assume that this 
time t FR is much shorter than the experimental time t, 
so that many loops are generated by one segment before 


Solid 4 He with finite density of superclimbing disloca¬ 
tions in a contact with superfluid reservoir is found to 
be, in general, characterized by bistability with respect 
to the syringe fraction. This feature is due to the in¬ 
terplay between three contributions: i) the chemical po¬ 
tential energy gain due to a transfer of atoms between 
two phases - solid and superfluid; ii) the energy of the 
deformation of dislocations needed to accommodate the 
transfer; and iii) the collective elastic energy. The con¬ 
trol parameter of the system are chemical potential and 
the density of the dislocations. At low densities of dis¬ 
locations the two fractions are very different from each 
other and, therefore, the transition between them can be 
viewed as strongly 1 st order transition with significant 
hysteresis. There is a similarity between this and liquid- 
gas transitions, with the exception of no critical point in 
the first case. 

It is highly likely that the syringe effect observed in 
Refs.dHJ is essentially in the overcritical regime. In this 
regime the equilibrium syringe fraction is given by the 
liquid type isochoric compressibility despite that the lin¬ 
earized response may predict much smaller values. 

The are two major scenarios for the instability. First, 
the vycor-solid boundary can be a source of the super¬ 
climbing dislocation loops entering the bulk. Its dynam¬ 
ics is characterized by the ballistic and dissipative stages 
which can be ohmic or strongly non-linear in the flow ve¬ 
locity. Each regime is characterized by specific powers of 
the bias and time at the initial stages of the evolution, 


Eqs.( |37|52|53 ). 
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Second, there is also an option for the bulk syringe ef¬ 
fect where the accumulated fraction is distributed evenly 
through out the bulk. In its turn, the bulk scenario 
can proceed in two ways - through the Bardeen-Herring 
generation of the prismatic loops or through the heli¬ 
cal instability of screw dislocations. The accumulated 
fraction in the case of the Bardeen-Herring instability 
is determined by constant rate dependencies (58), (60) 


in the non-linear and ohmic regimes, respectively. The 
non-linear regime (581 turns out to be showing the same 


type of the dependence of the syringe rate on the bias as 
the flow rate through the sample observed in the Ref. 1 . 
In this regard, it is worth mentioning a possibility that 
the flow through solid may not actually be taking place 
through a static network of dislocations percolating be¬ 
tween both vycor electrodes. Instead, the loops gener¬ 
ated during the Bardeen-Herring cycles may eventually 
be moving between two vycor electrodes. These loops 
are mobile due to the superflow along their rim and can 
serve as ’’vehicles” transporting the mater across a sam¬ 
ple. Center of mass speed V cm of such a loop is locked 
to the speed of the superflow V along its rim by a sim¬ 
ple geometrical relation V cm ~ Vb/R stemming from the 
matter conservation. Experimental studies of the actual 
bias-time-temperature dependencies of the syringe frac¬ 
tion and rates are needed to see if any of the above sce¬ 
narios take place. 

One of the key questions to answer is about the nature 
of the T-dependence observed in the non-linear regime of 
the flow rate 1 . Similar dependence w as also observed in 
a different setufP. Eqs.( 52|54|58|60 ) contain the super¬ 
fluid density p s in the corresponding powers as overall 
factors. To what extent the observed temperature de¬ 
pendence can be attributed to these factors remains to 


be seen. One possibility could be that the superfluid¬ 
ity along the cores is strongly affected by structural ex¬ 
citations of dislocation - kink^® and jogs - so that as 
T increases these excitations suppress the overall super¬ 
fluid density p s in the cores and, thus, reduce the total 
flow rate. It should also be mentioned that superclimb¬ 
ing dislocation does not fit exactly into the paradigm of 
Luttinger liquid because its excitation spectrum is not 
linear in the momentum^. To what extent this feature 
may modify the results ( 46|48 ) is an open question too. 

The ’’smoking gun” evidence for truly superfluid flow 
would be the detection of the ballistic jump in the syringe 
fraction which is a decreasing function of the bias. Some 
jumps have been observed in Ref. 2 3 4 5 6 7 . Thus, their detailed 
study is of crucial importance. 

The main assumption of this work is that density of 
superclimbing dislocations is low and a sample size is 
large enough so that the equilibrium for the generated 
loops is achieved at typical sizes smaller than sample size. 
If this condition is not satisfied, as it could be the case for 
very small sampled, a completely different scenario may 
take place: the conducting network may be created by 
dislocations proliferating directly between the electrodes. 
In this case the actual dynamics may be controlled by 
changing number of the conducting pathways, that is, 
balanced by the pathways creation and exiting from a 
sample. 
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